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ABSTRACT 

We report the discovery of a photometric variability in the bright T2.5 brown dwarf SIMP 
J013656. 57+093347. 3. Continuous 7-band photometry has been obtained for several hours on four differ- 
ent nights. The light curves show a periodic modulation with a period of ~2.4 hours, a peak-to-peak amplitude 
of ~ 50 mmag and significant night-to-night evolution. We suggest that the light curve modulation is due to 
the brown dwarf's rotation and that the longer term variations come from surface features evolution and/or 
differential rotation. We obtained complementary observations over a single night in the J and bands; the 
object displays correlated photometric variability in both bands, albeit with smaller /iTs-band amplitude. The 
ratio of the /Tj and J variability amplitudes puts strong constraints on the physical mechanisms at play. Based 
on theoretical models of brown dwarf atmospheres, our results suggest that the atmosphere of SIMP0136 is 
comprised of both grain-free and colder (by ^ 100 K) grain-bearing cloudy regions. This discovery, and its 
interpretation, provide a natural explanation of the so-called /-band brightening. 
Subject headings: Stars: individual(SIMP J013656. 5+093347) - stars: low-mass, brown dwarfs 



1. INTRODUCTION 

The study of brown dwarfs (BD) is a very active field in 
astronomy. Since the discovery of the first ones less than 15 
years ago (Rebolo et al. 1995; Nakajima et al. 1995), much 
has been learned on the nature of those objects. While spec- 
tacular advances in modeling their atmospheres have been 
made, understanding the physical properties of BDs strad- 
dling the L/T transition remains one of the most challenging 
problems. Among the unresolved issues is the intriguing in- 
crease of the 7-band luminosity with decreasing temperature 
over the 1400-1100 K interval, just following the transition 
from L to T type. This odd behavior, the so-called /-band 
brightening (Vrba et al. 2004), is surprising as one would 
naively expect near-infrared brightness to decrease with tem- 
perature. The brightening coincides with the disappearance 
of silicate grains from the brown dwarfs' atmospheres and is 
most likely a consequence of this disappearance, but a self 
consistent explanation has yet to be found. While dusty mod- 
els containing no grain sedimentation best reproduce the be- 
havior of late-M to mid-L dwarfs, condensed (clear) models, 
where the grains are entirely removed from the atmospheres, 
accurately describe mid-to-late-T dwarfs. Models including 
grain settling and precipitation (Allard et al. 2003; Burrows 
et al. 2006) have been able to reproduce a flattening of the 
7-band magnitude over the problematic temperature interval, 
but no single set of physical parameters can entirely reproduce 
the observed early-T dwarf sequence. The fundamental as- 
sumption of these models, that a single set of physical param- 
eters can describe the photosphere of T dwarfs, might not be 
accurate. The problem is further complicated by the possibil- 
ity that, to some extent, the /-band brightening may be caused 
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by a high frequency of unresolved binaries over this spectral 
type interval (Liu et al. 2006). Tsuji (2003, 2005) proposed 
a set of models with dust settling that describes the whole 
L and T sequence with a parameterization of the cloud deck 
thickness. They introduced a critical temperature, T^r, below 
which the dust no longer remains in the atmosphere in the 
form of clouds and is removed from the atmosphere through 
gravitational settling. They suggest that the /-band brighten- 
ing may be explained by the fact that the L-T sequence is not 
necessarily a Teff one but that it must instead be understood as 
arising from variations of both Teff and Tcr. Aside from this 
possibility, it has been suggested that clearings in an other- 
wise cloudy atmosphere may be responsible for the observed 
7-band brightening (Burgasseret al. 2002). In such a scenario, 
where the atmosphere is not uniform, photometric variability 
could be observed due to the rotation and/or evolution of the 
cloud pattern. 

There have been many observational searches for variabil- 
ity in ultracool dwarfs (M7 and later) in the last few years. 
Early searches concentrated mostly on L dwarfs in the / band 
and clearly established that roughly half of these objects dis- 
play variability on timescales ranging from tens of minutes 
to weeks, with amplitudes of a few tens of mmag (Tinney & 
Tofley 1999; GeHno et al. 2002; Bailer- Jones & Mundt 2001; 
Terndrup et al. 1999; Martin et al. 2001). Several studies 
suggest that variability also extends to BDs at the L/T tran- 
sition. Enoch et al. (2003) reported relatively large vari- 
abihty {AK^ 190-420 mmag with uncertainties of 40-50 
mmag) for 3 transition objects (L7, Tl and T5) with some ev- 
idence of periodic (1 -3 h) modulation. Other groups (Clarke 
et al. 2008; Koen et al. 2004, 2005) have detected variability 
for a few objects at a level of ~ 15 mmag. Recently, Goldman 
et al. (2008b) described a marginal detection of variability 
for the T2 dwarf SDSS 1254 using near-infrared spectroscopy. 
Finally, Morales-Calderon et al. (2006) tentatively detected 
photometric modulations (^^10 mmag at 4.5/im) in the ther- 
mal infrared for two late-L dwarfs. Overall, while variability 
seems to be established for L/T BDs, the significance of these 
detections remains marginal and no clear correlation between 
photometric bands has been reported. 
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TABLE 1 

Properties of the reference stars. 



2MASS designation 






Sep. from 
SIMP0136 


Comment 


JO 1365662+0933473 


13.45 


12.56 




SIMP0136 


JO 13707 17+0935071 


13.43 


12.54 


3.5' 


Comparison star 


J01364524+0931540 


11.82 


11.54 


3.7' 


Reference star 


JO 1370290+0937 184 


12.72 


12.35 


4.3' 


Reference star 


J01371365+0933327 


11.85 


11.18 


4.7' 


Reference star 



»2MASS(Skrutskie et al. 2006) point source catalog. 



The T2.5 brown dwarf SIMP J013656.57+093347.3 (Arti- 
gau et al. (2006), SIMP0136 hereafter), the second brightest 
known T dwarf after e Indi Bab, is an ideal target for study- 
ing the intriguing physical properties of BDs straddling the 
L/T transition. It falls in the middle of the 7-band brighten- 
ing and it does not appear to be a binary, down to separations 
of 0.2" and AK = 2 as inferred from Lick AO observations 
(Goldman et al. 2008a). In this paper, we report on the de- 
tection of significant near-infrared photometric variability for 
this object. The analysis of these observations not only pro- 
vides compelling evidence for the presence of clouds at the 
surface of this BD but also offers an attractive explanation for 
the 7-band brightening observed in early T dwarfs. 

2. OBSERVATIONS AND DATA REDUCTION 

The y-band observations were carried out at the Observa- 
toire du Mont-Megantic (OMM) 1.6 m telescope with the 
CPAPIR near-infrared camera (Artigau et al. 2004), featuring 
a 0.89" pixel scale and a 30' x 30' field of view. The obser- 
vations were obtained on 2008 September 16, 18, 19 and 21, 
under photometric conditions. As the CPAPIR point-spread 
function (PSF) would be undersampled under better than 1.8" 
seeing, and given that SIMP0I36 is a relatively bright tar- 
get, the camera was purposely defocused to produce 6-pixel 
FWHM PSFs. The individual exposure times were 14.8 s with 
a single coadd. The data were obtained without dithering and 
using a guide probe. To further minimize drifts of the field 
during the sequence due to flexure between CPAPIR and the 
guide probe in the telescope bonnette, a script was run in real 
time to send slow (~ 10 minutes), closed-loop corrections to 
the guide probe. The results from the first set of observations 
prompted us to take consecutive J and /Q-band photometric 
observations on the night of 2008, December 13. The data 
were also obtained at the OMM but using the SIMON near- 
infrared spectro-imager with a 0.46" pixel scale and a 7.8' 
field of view. The J and K^, sequences lasted 3.0 and 2.0 hours 
respectively with individual exposure times of 29 s in both 
bands. The observing strategy was the same as before. For all 
datasets, a high signal-to-noise (>10* e~/pixel) dome flat was 
obtained at the end of the night. 

All datasets were reduced using the same method. First, 
the images were dark subtracted and divided by the median- 
combined dome flat. Then boxes of 50 x 50 pixels centered on 
SIMP0136, on three nearby reference stars slightly brighter 
than SIMP0136 and a on comparison star closely matching 
the / and Ks of SIMP0136 (see Table 1). The same reference 
and comparison stars were used for both the CPAPIR and SI- 
MON observations. A handful of bad pixels, none of them 
falling on the PSFs, were present in the boxes; they were sim- 
ply interpolated. Aperture photometry was performed using 
radii of 2.5 times the 50% encircled energy radii of the PSF 
(typicaUy 5-6 pixels). The sky level was estimated by using 
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JD + 2 454 731 

Fig. 1. — Raw light curves for SIMP0136 and the comparison star on 
September 21, 2008. The modulation is clearly seen in SIMP0136 light 
curve. The slow modulation in the reference star light curve is seen for all 
objects and is due to the airmass change during the night (~1.95 at the be- 
ginning and end of the sequence and a minimum of 1.24 at transit). 
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Fig. 2. — The Ught curves of SIMP0136 for each observing night are shown 
on the upper panel. The lower panel shows the light curves of the comparison 
star. The light curves are binned over 5-minute time bins. 

the median value of all pixels beyond 1.5 times the photomet- 
ric aperture radius in each box. 

The above steps produced raw light curves for the target 
and each of the three reference stars. The raw light curves of 
SIMP0136 display a strong modulation unUke anything seen 
in the raw Ught curves of the reference or comparison stars 
(Figure 1). The raw curves of SIMP0136 and the comparison 
star were then corrected for systematic effects in the follow- 
ing manner. The light curve of SIMP0136 and the comparison 
star (shown in Figure 2) were obtained by dividing them by 
the mean normalized raw light curves of the 3 reference stars, 
see Figure 2. We obtained between 460 and 867 individual 
images per night and the photometric points shown in Figure 2 
are binned over 5-minute intervals (11 to 14 individual mea- 
surements). Uncertainties per individual image ranged from 
14 to 19 mmag for both SIMP0136 and the comparison star, 
corresponding to a '--^ 5 mmag uncertainty over a 5-minute in- 
terval for both objects. A similar analysis was done for the 
December 13 data and the results are shown in Figure 3. 

3. DISCUSSION 

The photometric variability of SIMP0136 is clearly visible 
from Figure 2. The variation reached a peak-to-peak ampli- 
tude of AJ ^ 50 mmag. The December data (see Figure 3) 
clearly show that the J and light curves are correlated 
with one another, the variations being more pronounced in 
the former with an amplitude ratio of ^f- = 0.48 ± 0.06. The 
September 21 light curve, covering 7.5 h, shows a periodic, 
double-peaked signal. A very similar modulation is also ob- 
served on September 19. As a rough estimate of the period 
can be readily obtained from the September 21 light curve, 
and as there is an important second harmonic component in 
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Fig. 3. — J and Ks band light curves obtained on December 13, 2008. 
The upper panel shows the complete light curves of SIMP0136 and the com- 
parison star while the lower one shows the light curves of SIMP0136 folded 
to a 2.895 h period. The modulation observed in September 2008 is still 
clearly present and observed in both bands. The amplitude of the mod- 
ulation is correlated with the /-band one but with an amplitude 0.48 it 0.06 
times smaller The solid line in the J-band phase-folded light curve shows 
the smoothed curve and the continuous line in the A^s-band phase-folded light 
curve is the same line multiplied by 0.48. The phase-folding relative to the 
September 2008 observations is arbitrary. 

the data set, we determined the precise period by fitting a 
two-harmonic signal using a Levenberg-Marquardt algorithm 
rather than using a Lomb-Scargle period-search algorithm 
(Lomb 1976; Scargle 1982), which assumes a sinusoidal sig- 
nal and is better suited for sparsely sampled datasets. We 
derived a 2.3895 ± 0.0005 h period using either the Septem- 
ber 19 and 21 datasets. The evolution of the Ught curve of 
SIMP0136 prevents us from deriving a more accurate period 
as the light curve cannot be folded over several nights. Be- 
sides, differential rotation probably renders pointless the at- 
tempt of estabhshing a much more precise period. 

A 2.3895 h period would correspond to a vsinj of up to 
^ 45 km/s depending on the object's rotation axis inclination 
with our line of sight, consistent with the fastest rotators in 
the Basri et al. (2000) and Zapatero Osorio et al. (2006) sam- 
ples. Only a handful of brown dwarfs have a similar rotation 
period. The L2 dwarf Kelu- 1 has a rotation period of less than 
2 hours (Clarke et al. 2002). The young M8.5 dwarf S Ori 45 
displays a photometric (in the J and / band) modulation with 
a ~ 2.6 hours period but this periodicity was not detected in 
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Fig. 4. — Phase-folded light curve of SIMP0136 for each observing night 
assuming a period of 2.3895 h. The continuous line in each epoch's graph 
represents the smoothed light curve observed on 2008 September 21. An evo- 
lution of the light curve is clearly observed at successive epochs. The double- 
peaked feature clearly detected on September 19 and 21, and marginally on 
September 16, has disappeared on September 18. 

later observations despite hints of a yet faster (30-45 minutes) 
modulation (Zapatero Osorio et al. 2003). Strong discrete pe- 
riodograms peaks have been reported for 2M1 145 at periods 
ranging from 2 to 3.5 hours (Bailer- Jones & Mundt 2001) but 
the light curve of this object cannot unambiguously be folded 
to a single period. 

The clear evolution of the light curves from night to night 
in Figure 4 could be due to time- varying clouds or storms in 
the atmosphere of SIMP0136. In that case, observations over 
several days can give a time scale for the evolution of those 
clouds/storms. The evolution of the modulation could also 
be attributed to a differential rotation. If such an assumption 
is made, a significant phase difference between the equator 
and the higher latitudes has to be present within our five-day 
baseline (here we assume that the December observation is 
entirely decorrelated). Assuming a (maximum) phase differ- 
ence of 180° between the equatorial and high-latitude regions 
of S1MP0136 over the 5-day-long observing sequence would 
imply a 1% differential rotation. This value is comparable 
to Jupiter showing a ~ 0.8% differential rotation per period 
between its equator and high latitudes. It is thus possible that 
the evolution of the light curves in Figure 4 is the result of dif- 
ferential rotation. This hypothesis will require a much better 
temporal sampling to be confirmed. 

The relative variations at / and can be used for constrain- 
ing the physical processes at play. Using the AUard et al. 
(2003) and Burrows et al. (2006) models, assuming an effec- 
tive temperature of 1200 K, logg = 5.0 and solar metallicity, 
we derived for three different variability scenarios: 1) 
dusty clouds on a clear atmosphere with both phases at the 
same effective temperature, 2) inhomogeneous surface tem- 
perature, for either a dusty or a clear atmosphere and 3) dusty 
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clouds on a clear atmosphere with both phases at different 
temperatures. Here dusty and clear refer to the settled and 
condensed models of Allard et al. (2003) and to the cloudy 
and clear models of Burrows et al. (2006). For scenarios 2 
and 3 we further considered the models of Tsuji (2003, 2005), 
the case where dust clouds are completely absent from the at- 
mosphere is presented only as a limiting case (called case C 
in Tsuji (2003, 2005)), it assumes a Tcr equal to the conden- 
sation temperature of the dust grains. 

For the first scenario, both sets of models predict that, for an 
atmosphere of constant temperature, the J and band light 
curves should be anti-correlated (i.e. ^^<0). This can easily 
be understood as the cloud-free models are brighter than the 
cloud-bearing models at J and inversely at K^. Thus if an 
object evolves from, say, a 50% cloud-free and 50% cloud- 
bearing atmosphere toward a 51% cloud-free and 49% cloud- 
bearing atmosphere, it will brighten in / and dim in K^. 

In the second scenario, the variability is caused by surface 
temperature inhomogeneities in an entirely clear or dusty at- 
mosphere resulting in a time-variation of the disk-averaged 
temperature. Numerically, this is done by taking the ratio 
of the relative flux increases in J and Kg, for a small (i.e. 
Ar <C T) temperature change : 



A/ 



V 9T ) 



Fj 



m fk. 



(t) 



Using all models, we derived a ^f- ranging from 1.2 to 3.4, 
the lower and upper limits being set respectively by the cloudy 
and settled models of Allard et al. (2003). The Tsuji (2003) 



models using a Tc, = 1800, yield a 



Ay 



'2.0. 



The observed = 0.48 ± 0.06 allows us to exclude the 
first two variability scenarios. We now consider the third case 
whereby clouds have a different temperature from the sur- 
rounding cloud-free atmosphere. Let / be the cloud filling 
factor and A/ the variation in cloud coverage. We assume 
that clear (cloud-free) regions are at a temperature T = 1200 
K, appropriate for a T2.5 dwarf, and that clouds have a tem- 
perature of r + Ar. Let F^IT] and F^T] be the fluxes in the 
photometric bandpass m as a function of temperature for the 
dusty and clear models. The variability is caused by the vari- 
ation in cloud coverage (A /) as the brown dwarf rotates. It 
can be shown that the resulting variability Am in a photomet- 
ric bandpass is 



Am = 2.5 log 



FdT + AT]f+F,[T]{l-f) 

Fd[r+Arj(/+A/)+FeLrj(i-/-A/) 



The ratios for all models and different values of / 
and Ar are shown in Figure 5. As shown in this figure, 
the resulting varies only slightly as a function of the 
assumed filling factor for clouds versus cloud-free patches. 
Also, the overall behavior between the Allard et al. (2003) 
and Burrows et al. (2006) models is consistent. The observed 
= 0.48 ± 0.06 can be reproduced by clouds between 80 to 
1 10 K colder than the cloud-free atmosphere. The results are 
mildly dependent on the temperature of the clear regions; a 
r=1300 K would increase the derived AT to values ranging 
from 125 to 150 K. This variability scenario also gives a ^ 
ratio allowing for a direct conversion of the light curve ampU- 
tude (AJ) into an estimate of the variation of the cloud filling 
factor (A/). The derived ^ range from 1.5 to 2.3 depend- 
ing on the assumed filling factor / (between 0.25 and 0.75). 
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Fig. 5. — Ks versus J band variability as a function of the temperature 
difference between the cloudy and cloud-free atmosphere patches. The hori- 
zontal line corresponds to the observed and the thickness of the line to 
the imcertainties on that measurement. The Allard et al. (2003) models do 
not extend below Ar = -100 K as the settled models have only been com- 
puted down to 1 100 K. All models assume a temperature of 1200 K for the 
cloud-free atmosphere, logg = 5.0 and solar metallicity. 

This implies that for a 50 mmag variability amplitude in /, 
the patches involved in SIMP0136 variability cover between 
7% and 16% of the visible disk. This does not constrain the 
overall filling factor of the dusty versus clear regions, only 
time-differences. 

In the Tsuji (2003, 2005) models, the case C (equivalent to 
the clear models) is not expected to cover a large fraction of 
the surface, we therefore examined only the case / = 0.75 case 
where 75% of the surface would be covered by the Tcr = 1800 
model (equivalent to the dusty models). The observed 
implies a ~ lOOK temperature difference between the case C 
and cloud-bearing regions, consistent with the predictions of 
the Allard et al. (2003) and Burrows et al. (2006) models. 

Beside atmospheric phenomena, periodic modulation at a 
few percent level could easily be introduced by the presence 
of a transiting companion wifli a radius of 2-3 or by graz- 
ing transits in a roughly equal-sized binary. However, the 
observed night-to-night evolution of the light curves strongly 
disfavors any explanation of the variability in terms of tran- 
sits/grazing transits. The expected light curves would be very 
stable from period to period and evolve extremely slowly 
through the evolution of the orbital elements. Also, any transit 
scenario would imply an out-of-transit flat-top light-curve that 
is not observed. Furthermore, transit light curves show little 
wavelength dependence (apart from limb darkening effects) 
and one would expect a ~ 1 regardless of the geometry 
of the system. 

4. CONCLUDING REMARKS 

The observed variability of SIMP0136 suggests that atmo- 
sphere models attempting to describe its physical properties 
using a single set of parameters (dust settling rate, tempera- 
ture, gravity, etc) are bound to be incomplete as at least two 
different types of regions must be present on its visible disk to 
produce the observed variability. The idea that early-T dwarfs 
could be described as objects with more than one set of phys- 
ical properties and arbitrary filling factors is not new (Bur- 
gasser et al. 2002), but the detected variability and its inter- 
pretation in terms of colder clouds in a warmer atmosphere 
provides empirical physical constraints of these cloudy and 
cloud-free regions. 

It has been proposed previously that the 7-band brighten- 
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ing may result mainly from a decrease in the fractional cover- 
age of dust clouds against warmer, clear regions of the brown 
dwarfs' atmospheres as they progress from spectral type L to 
T. Indeed, a brightening in 7 is a natural outcome of this as 
warmer clear regions are brighter than cooler dusty regions 
in /. Our observations not only bring empirical evidence for 
this scenario, but also provide an estimate of the temperature 
difference between the cloudy and clear regions of the atmo- 
sphere. The difference of ~100 K suggested by our results 
would lead to a brightening of 0.64 mag in J from a com- 
pletely cloudy to a completely clear atmosphere, assuming no 
temperature change for each respective model; this amplitude 
is similar to that of the observed /-band brightening. In real- 
ity, the L/T transition is likely accompanied by a mild change 
in the temperature of both types of atmosphere - cloudy and 
clear - but the overall behavior would be similar. The exact 
shape of the /-band brightening could potentially be repro- 
duced by adjusting, as a function of spectral type, both the 
temperature of the clear regions and the difference of temper- 
ature between the clear and cloudy patches. 

Considering that most searches for variability in early-T 
dwarfs resulted in non-detections or marginal detections de- 
spite sensitivities sufficient to readily detect variability at the 
level observed for S1MP0136, the following question natu- 
rally arises : is SIMP0136 an oddball object with unusual 
physical properties ? or is it representative of early-Ts ? The 
non-detection of variability in other early-Ts does not dis- 
prove the cloudy versus clear hypothesis: cloud structure on 
these objects could be either too small to produce significant 
photometric variability or distributed in bands without longi- 
tudinal inhomogeneities. The near-infrared photometric prop- 
erties of SIMP0136 are typical of early-T dwarfs (See Figure 
2 in Artigau et al. (2006)) and its proper motion is consis- 
tent with that of a thin disk object; it is therefore not expected 
to have unusually low metaUicity and/or high surface gravity. 
The only slightly unusual property of SIMP0136, besides its 
photometric variability, is its relatively rapid rotation. Deter- 
mining whether rapidly rotating brown dwarfs are more likely 
to feature larger-scale atmospheric structures would require 
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global circulation models; this could be verified observation- 
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